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Abstract 

This paper proposes a bottom-up sequence of modeling steps leading to a nanoscop- 
ically informed continuum, and as such macroscopic, theory of single-walled car- 
bon nanotubes (SWCNTs). We provide a description of the geometry of the two 
most representative types of SWCNTs, armchair (A-) and zigzag (Z-), of their 
modules and of their elementary bond units. We believe ours to be the simplest 
shell theory that accounts accurately for the linearly elastic response of both 
A- and Z- CNTs. In fact, our theory can be shown to fit SWCNTs of whatever 
chirality; its main novel feature is perhaps the proposition of chirality-dependent 
concepts of eff'ective thickness and effective radius. 
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1 Introduction 

Carbon has a large variety of allotropic forms, natural and artificial: graphite, 
diamond, fuUerene, carbon fibers, carbon nanotubes (CNTs), and others. The 
relative C-C bond complexes, diverse as they are by spatial dimensions and 
geometry, correspond to diverse hybrid electron states, made possible by com- 
bination of the electrons occupying the most external shells in a C atom; the 
nature and the strength of its C-C bond complex decide the mechanical prop- 
erties of a given allotrope at any scale. 
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This paper proposes a bottom-up sequence of modeling steps leading to a 
nanoscopically informed continuum, and as such macroscopic, theory of single- 
walled carbon nanotubes (SWCNTs). Recently, modeling procedures of the 
same type have been adopted by various authors [13 IH HH UHl HOI M 121] , with 
application to CNTs and polymer composites in mind. Related multiscale mod- 
eling papers, such as [T] and (TB], cover a variety of tools for studying both quan- 
tum and classical models of atomic systems in crystalline phase, and propose 
techniques to connect the microscopic models to the continuum limits. They 
suggest that, for non-periodic systems, the key challenge lies in carefully defin- 
ing average energy at the microscopic level and connecting it to a macroscopic 
energy. Another related paper is [TH], that studies the mechanics of materi- 
als with periodic micro-structures; here, lattice structures are divided into two 
classes: truss materials modeled as beams with axial and bending stiffness, and 
frame materials modeled with pin-joints at the nodes. 

There are of course differences in concept ingredients and developments 
in the cited papers, but their goal is one and the same: to deduce the tar- 
geted macroscopic mechanical behavior from the available nanoscopic chemical- 
physical information. In our present paper, a bridge between the nanoscopic 
and the macroscopic viewpoints is created by means of an intermediate pas- 
sage, that we term meta-nanoscopic, in a fashion that we now describe. At 
the nanoscopic scale, the bonding and non-bonding energies keeping CNTs to- 
gether are evaluated in the fashion of molecular mechanics (MM) (Section [5]). 
In the following section, a mechanical caricature of a SWCNT as an orderly 
arrangement of pin-jointed sticks and (axial and spiral) springs is drawn, within 
the framework of discrete structure mechanics (DSM). This requires a prelim- 
inary careful description of the geometry of the two most representative types 
of SWCNTs, armchair (A-) and zigzag (Z-), of their modules and of their el- 
ementary bond units; in such description, the nanoscopic parameters are the 
chirality index and the length of the C-C bond. The viewpoints of MM and 
DSM are connected by equating the energies per bond, for each of which an 
approximate quadratic expression is accepted: two types of bonding energies 
- hence two purely constitutive parameters, the stiffness constants of the two 
types of springs - enter the expression for the DSM energy of a bond unit; they 
are taken equal to the two corresponding nanoscopic stiffness parameters, re- 
garded as a known input from MM. This simple measure opens the way to the 
posing and solving, within the framework of DSM, of a number of equilibrium 
problems for SWCNTs (Section El). 

The step from discrete to continuous structure mechanics (CSM) is harder. 
To take it, a DSM model of SWCNT has to be assimilated to a suitable model 
taken from the library of CSM. The main difficulty in doing so resides in the 
number and nature of the parameters needed to specify a CSM model: no matter 
what rod or shell theory one chooses, the relative stiffness notion consists of a list 
of more than two combinations of both constitutive and geometric parameters, 
the latter reflecting theory-specific concepts of thinness and slenderness. Now, 
while there is a natural nanoscopic notion of slenderness for SWCNTs, to assess 
their thinness is a controversial issue, because a preliminary estimate of their 
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effective thickness is required (these matters are dealt with in Section HTT)) : and, 
no doubt both slenderness and thinness of a SWCNT depend on its chirality. 
Therefore, one is faced with the problem of determining the chirality-dependent 
parameters of the continuum theory of his choice in terms of only two nanoscopic 
constitutive parameters and in the absence of an unequivocal nanoscopic concept 
of thinness. 

For reasons that are discussed at length in [8] and recapitulated in Section [5l 
the CSM theory we adopt is a linearly elastic theory of orthotropic shells, where 
the independent constitutive parameters are four and the geometric parameters 
two. To determine all those parameters, firstly we solve the system of equilib- 
rium equations ruling the axial-traction and radial- loading DSM problems for 
both A- and Z- CNTs, and the torsion problem for ACNTs ( Sections l4Jl - |43)) : 
interestingly, the solutions we find hold for whatever slenderness. Secondly, 
we import from |5] the analytic solutions there found for the same equilibrium 
problems, when set within the framework of the above mentioned shell theory. 
While the torsion solution is insensitive to slenderness, the axial-traction and 
radial-loading solutions are; for simplicity, we confine attention to the case of 
slender shells. Thirdly, for each of these boundary-value problems, both in the 
discrete and in the continuous formulation, we regard the involved nanotubes as 
cylindrical probes, whose relevant stiffnesses and Poisson-like contraction ratios 
we define in the fashion of standard rod theory (Sections 16.11 and 16.21 respec- 
tively). Finally, in Section [6. 3[ the sought-for parameters are found by equating 
the continuous and discrete stiffnesses and contraction ratios; since we have an 
unusually rich crop of solutions at our disposal, this can be done in more than 
one way, with convincingly consistent results. 

We believe ours to be the simplest shell theory that accounts accurately for 
the elastic response of both A- and Z- CNTs. In fact, modulo a modest gener- 
alization to be presented elsewhere, our theory can be shown to fit SWCNTs of 
whatever chirality; its main novel feature is perhaps the proposition of chirality- 
dependent concepts of effective thickness and effective radius; for a discussion 
of this issue, and others, we refer the reader to our conclusion section. 

2 Nanoscale Energies 

At nanometer scale, atomic interactions are modeled by either quantum mechan- 
ics (QM) or molecular mechanics (MM), in both cases by striving to capture how 
and to what extent the system's energy varies with changes in atomic positions. 
QM can describe rigorously the electronic structure of a material complex, but 
its computational cost quickly becomes prohibitive as the number of atoms in- 
volved increases. MM is based on the Born-Oppenheimer approximation for 
the Hamiltonian of a collection of heavy and light particles, a cornerstone of 
computer simulations. Consistent with this simplifying assumption, the total 
energy U of such a material complex is expressed as a sum of two terms: 
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where Ub and Unb, denote, respectively, the bonding energies and the non- 
bonding energies. The latter account for van der Waals and Coulomb interac- 
tions. However, when modeling SWCNT, the contribution of Coulomb interac- 
tions is usually neglected since interactions are predominantly between essen- 
tially neutral carbon atoms. Within the same modeling context, the former 
consist of four parts: 

Ub^Up + Ue + U^ + Ur, 

where Up,Ue, Uuj and Ur denote the energies respectively associated with stretch- 
ing (of a covalent bond), angle variation (between two covalent bonds), torsion 
(around bonds), and the so-called improper torsion (see e.g. 19 J. The third 
and fourth contributions are considered negligible when compared with Up and 
Ue, whose harmonic approximations are: 

Up = lJ2''p(p~p"'')'^ 

C/, = i5]fc,(^?-0'- 

Thus, the only bond-stiffness constants of importance are kp and kg; they can 
be obtained by ab initio QM evaluations or fitted to experiments. The values 
we shall employ in our numerical computations, 

kp = 742 nN/nm and kg = 1.142 nN x nm, 

are the same as in [4] and |24) . 

3 The Geometry of SWCNTs, from A to Z 

CNTs were discovered in 1991 by S. lijima |12j. They look like right cylinders 
with approximately emispherical caps due to the actual manifacturing proce- 
dures, and can be single- or multi-wall. In imagination, SWCNTs can be ob- 
tained by rolling up into a cylindrical shape a graphene, that is, a monolayer 
flat sheet of graphite that can be visualized as a two-dimensional lattice with 
hexagonal unit cell. There are many ways to roll a graphene up, sorted by 
introducing a geometrical object, the chiral vector: 

X = nai + ma2 

(Fig. [T]). Having fixed ai and a2, the pairs of integers {n,m) specifies the 
chirality of a CNT; (n, 0)— and (n, n)— nanotubes are called, respectively, zigzag 
and armchair (Fig. [SJ roll- up axes are chosen perpendicular to chiral axes)lf| 

Chirality influences all geometrical and mechanical properties of a CNT; for 
example, the radius of the cylinder on which a (n,m)— SWNT sits is: 

Po =?o("-)™) ■= 77- V3(n2 +nm-\-m'^)s, (1) 

ZTT 

•^The chiral angle is defined to be x(n,m) = arctan(\/3 ^^2n )■ thus, in particular, the 
chiral angle of an armchair CNT is 7r/6 radians. 
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(n,0) zigzag 



Figure 1: Graphene and its zigzag and armchair chiral axes. 



where s the side length of a graphene ceh; and, chirahty influences the effective 
extension and shear moduli and the buckling strain [51 131 [3 [531 ISI] ; needless 
to say, for multi-wall carbon nanotubes (MWCNTs), the chirality of the inner 
waUs is hard to inspect. 

In this paper, we concentrate on SWCNTs, armchair and zigzag; for short, 
we designate them by ACNTs and ZCNTs, respectively. The geometry of 
their undeformed reference configurations is described just below. To get ac- 
quainted with their size, it is convenient to start from a graphene cell: its 
side length s equals the C-C bond length, about 0.142 nm (that is, roughly, 
the diameter of a C atom); consequently, two opposite C atoms are spaced by 
(1 -I- 2cos7r/3)s — 0.283 nm = the cell's diameter; and, equally placed C atoms 
belonging to adjacent cells are spaced by \ai\ = \aQ\ = (2sin7r/3)s = 0.2456 
nm, the distance of two parallel sides of a cell (to fix ideas, the spacing between 
two adjacent cylinders in a MWNT is about 0.34 nm). 



Figure [31 is an atom-and-bond cartoon of the module of the ACNT having the 
smallest radius one can think of (n = 2), in its reference configuration. Whatever 
n, the geometrically necessary radius, that is, the radius of the cylinder on which 
the centers of the C atoms are placed, is: 



3.1 Module and bond unit of an ACNT 




n s 



(2) 



(cf. ()22|) ): the total number of atoms per module is: 



AM 



= 4 X 2n = 8n. 
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ARMCHAIR CNT 



ZIGZAG CNT 



chiral 
axis 






Figure 2: Graphenes to be rolled up into armchair (left) and zigzag (right) 
single-wall CNTs. Note the orthogonality of zigzag (red) and armchair (blue) 
atom sequences. 



distinguish two types of sticks, a horizontal and b oblique, and two types of 
angles between sticks, denoted by a and /?; in a module, there are: 



(3) 



Three adjacent sticks, two of type b and one of type a, form an armchair bond 
unit (ABU) . We make the geometry of an undeformed ABU precise with refer- 
ence to Figures m and [5l where atom A is sitting at the origin of an orthogonal 
cartesian frame of unit vectors c,;, with Ci = b. Note that 



a 


= 4 X 71 




type a sticks; 




= 4 X 2n = 


8n 


type b sticks; 




= 4 X 2n 


8n 


type a angles; 




= 4 X 4n = 


16n 


type /? angles. 



cos/3 



' cos — cos 7, 
2 ' 



7 



2n' 



(4) 



and that the three next neighbors of A sit at the ends of the relative position 
vectors: 

b = ci. 



s ^pc 



s ^Pd 



a 

c = cos — a - 
2 

a = cos — a 
2 



a 

sm-cg, 

. a 

- sm-C3, 



where the unit vector 



is orthogonal to C3; moreover. 



' COS7C1 + sin7C2 



- cos 27 Ci -t- sin 27 C2 . 



(5) 

(6) 
(7) 
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Figure 3: The module of a (2, 2)-CNT, rolled and not. 




Figure 4: The geometry of an ABU: front view. 



Remark. One may wonder why we have not set a = 27r/3 in the formulas 
just above. This is in fact the case for the undeformed configuration, when we 
deduce from ^ that 

/I TT \ 

P — (3 (n) = TT — arccos f - cos ^ j • (8) 

Our reason is that (U)-^ holds also when the applied loads induce deformations 
that, while changing angles of both types a and /3, do not change angle 7 and 
preserve the mirror symmetry with respect to the plane through A orthogonal 
to C3 of the position vectors of atoms C and D (e.g., these conditions are met 
in the axial-traction problem we are going to study in Section 14.11 but they are 
not in the torsion problem studied in Section 14. 3p . We record here for later use 
an expression of A/3 as a function of Aa that follows by differentiation of (jH)-^ 
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A B 
Figure 5: The geometry of an ABU: top-down view. 



while keeping 7 fixed: 



sni ^ ^ „ 

A/3 = — J (a, /3,7)Aa, 5 (a, /3, 7) := — COS7 = — — -C0S7. (9) 



73 



2sin/3 



4sin/3 



(cf. eq. (14) of 



3.2 Module and bond unit of a ZCNT 

The module of a ZCNT of smallest radius [n — 2) is drawn in Figure [SI Figure 
[7]is relative to the cases n = 4 and n = 10. The cylinder's radius and the total 





Figure 6: The module of a (2,0)-CNT. 
number of atoms per module are, respectively, 



p^{n,0) -. p^{n) ^ ns and N^''=4xn = 4n. 

27r 



In each module, there are: 



2 X n = 2n 
2 X 2n = 4n 



^ * a 



jyZM = 2 X 4n = 8n 
N^^'^ 2 X 2n = 4n 



type a sticks; 
type b sticks; 
type a angles; 
type /? angles. 



(10) 



(11) 
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Figure 7: A (4,0)- and a (10,0)-CNT. 




D .1 



B 




Figure 8: The geometry of a ZBU: front view. 



Figures [8] and |9] illustrate a zigzag bond unit (ZBU) and its geometry. We 
see that the three next neighbors of the typical atom A, labeled B, C and D, 
sit at the following positions: 



s ^pb ■= b = cos if Ci — sin ip C3 



s ^pc C3, 



(12) 



s ^pd '■= d = cos (fi a — simp C3, 



with the unit vector a as in ([6]); note that a = ip + tt/2. It is not difficult to 
show that angles a, /? and 7 are related as follows: 



whence 



cos f3 = cos a — sin a cos 7, 



~ B (n) — arccos ( cos — 



(13) 
(14) 
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Figure 9: The geometry of a ZBU: top-down view. 



On repeating the considerations proposed in Remark 1, and under the same 
circumstances, differentiation of (IT51) yields: 



A , r, A , r, sin2a , -\/3 , 

A/3 = 5 a,^,7 Aa, 5 /?, ^ := 1 + C0S7 = -— — 1 + C0S7 . 

smp 2smp 

(15) 

Remark. Here and henceforth, to hghten our notation, we write 13 in the place 

of /? (n) and /3 (n); no confusion should arise, because the type of nanotube 
we deal with is always clear from the context. 



4 Nanoscopic DSM. Axisymmetric Equilibrium 
Problems 

We view a meta-nanoscopic mechanical model of a SWCNT as consisting of pin- 
jointed rigid sticks and linearly elastic springs, of two types: (1) axial, sensitive 
to changes in distance of the two C atoms sitting at the ends of the coaxial 
stick, of stiffness ka, (2) spiral, sensitive to changes in angle of the two sticks 
they are attached to, of stiffness fc^. Our mechanical model of an ABU looks as 
in Figure [TUl (left): the zigzag case is shown in Figure [TT] (left). 

In this section we apply the principles and methods of DSM to solve explic- 
itly a number of axisymmetric equilibrium problems for SWCNTs, namely, axial 
traction and radial loading for both A- and Z- CNTs, and torsion for ACNTs. 
Problem by problem, different symmetries are presumed that simplify consid- 
erably the solution process: to guess those symmetries correctly is at times not 
trivial; once they are detected, it is a matter of elementary vector algebra. For 
this reason, in each of the subsections here below, we relegate the ancillary 
computations to the Appendix. 



4.1 Axial traction of an ACNT 

This problem has been dealt with in [H [151 El] , for zigzag CNTs as well. Figure 
[T2I is meant to give a helping hand to visualize an ABU and its immediately 
adjacent sticks when, on applying Euler's Cut Principle, that ABU is required to 
be an equilibrated part of a nanotube subject to axial traction and, accordingly, 
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Figure 10: A stick-and-spring caricature of an ABU. 

is regarded as acted upon by a balanced system of forces and couples, due to 
its interactions with the adjacent sticks. We assume that forces /b, fc, and fo, 
such that 

fB+fc+fD = 0, (16) 

are applied at points B, C , and D. In view of the peculiar symmetries of the 
axial traction problem, we presume that 

fB ■ 0, 

/c • C3 = -fo ■ C3, ^^^^ 
fc ■ Ca = fo ■ Ca = --fs ■ Ca (a = 1, 2) 

(note that, with these presumptions, p6)) is satisfied). 

In addition to forces, at both ends of each stick a number of couples are 
applied, which we distinguish into external and internal with reference to the 
Euler's cut which singled out the ABU. With a view toward writing moment 
balances, we stop and make our vectorial representation for couples unambigu- 
ous, with the help of Figure [TS] The figure alludes to a linearly elastic spiral 
spring of stiffness kg and rest angle ??, connecting two pin-jointed rigid sticks di- 
rected along the unit vectors u,v. When the angle changes by |At?|, the spring 
exerts a couple 

T=TUXv, T ^ ksA-d\u X v\^^ , 

on the stick along u and a couple T = — t it x t; on the other stick; thus, for 
A-d < 0, a clockwise (counterclockwise) couple is exerted on the stick along u 
{v). For example, upon loading the nanotube, the stick AC has external couples 
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Figure 11: A stick-and-spring caricature of a ZBU. 




Figure 12: An ABU and its adjacent sticks. 



Ta, applied in C and internal couples Tq,, applied in A, with 

Ta = fa (— c) X m, — fc^Ao! |c X m\~^ , 
T0 = -f0 e X (-c), T/3 = fcsA/3 |e x c|"^ 

and 

Ta = Ta C X d, Ta = kgAa |c X d\~^ , 

Tji = —Tp b X c, Tfj = ksA/3 \b X c\~^. 

With a further application of the Cut Principle, the moments of the force 
and couple system acting on stick AC are balanced with respect to the pole A 
if it so happens that 



(18) 



(19) 



PC X fc + Ta + T0 + Ta + T0 = 

where we presume that 

Aa = Aq, 'Ap = Aj3 



(20) 
(21) 
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Figure 13: The couple associated with a spiral spring. 



(the problem's symmetries exonerate us from laying down the relation for stick 
AD, which would not add any relevant information). Likewise, moments of 
force and couples system acting on stick AB are balanced with respect to A 
provided that 

Pb ^fB = 0, 

whence 

/b • C2 = 0. (22) 



Algebraic manipulations to be found in Section ISTT] lead us to conclude that the 
changes in length of sticks of type a and b are: 

^«=%-^-0' A6=:^=sin^//fc, = ^//fc„ (23) 
moreover, as to angle changes, we have: 



sinf / 1 /tanf \^\ ^ , / 3 , 



4tan;3 \ 2Vtan;3y / 8tan/3 V 2tan^/3^ 

(24) 

where angle f3 depends on the chirality index n as specified by ([8]) (recall Remark 



4.2 Axial traction of a ZCNT 

With reference to Figures [M] and [U let a system of forces satisfying ([TB]) be 
applied at points B, C, and Z? of a typical ZBU, singled out in Euler's manner 
to assess its equilibrium conditions as a part of a zigzag nanotube subject to 
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Figure 14: A ZBU and its adjacent sticks. 



axial traction. The built-in symmetries require that: 

fc = /C3, 

/b • C3 = /d • C3 = -i /c • C3 = -i /, (25) 

fo ■ a = fs ■ ci, /d • a X C3 = /b • C2. 

As to moment equilibrium, that of stick AC with respect to A is automatic, 
given the symmetries. Not so for stick AB: it must be that 

PB ^ fs + Ta + + Ta + Tf} ^ 0, (26) 

where 

Ta^Tab X C3, Ta = ksAa\b X C3 | 

Tq = fQ(-6) X (-C3), fa = ksAa\b X C3\ 
Tj3 = —Tp d X b, Tj3 — ksAl3\d x b\^^ 
T/3 = —ffj m X (—6), Tp = ksA(3\b x mp"^ 
with 6, d as in (IT2t and 



m — sina(cos7Ci + sin7C2) — cosacs. 
Equation (|26l) reads: 

b X (^sfs + 2tc,C3 + Tp{d - m)^ = 0, 
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where we presumed Aa — Aa and A/3 — A/3. After the algebraic manipulations 
confined in Section EUl we conclude that 



Aa = = //fca, A6= ^^-^ = -icosaj/fca = ^f/ka, 



and 



Aa = 



sm a 



4 1 + 2 



COS a sin^ a sin (3 



3 — cos /3 + cos 2q; (1 — 3 cos (3) 



8(5 + cos ;9) 



3sin/3 



4(5 + cos/3) 



(27) 



(28) 



4.3 Torsion of an ACNT 

This problem has been dealt with in [TSl [M] . The Euler cut in Figure [15] indi- 
viduates a double ABU, a five-stick complex that we require to be in equilibrium 
under a system of forces and couples applied at points C,D,E, and F. 




Figure 15: A double ABU and its adjacent sticks: front view. 

As to forces, on guessing the problem's symmetries, we assume that 

fc = htc + VC3, fo = -htc + VC3, ^^^^ 
fs = Me - VC3, fp = -Me - wca, 

with Ie the unit vectors orthogonal to C3 and tangent at points C and E, 
respectively, to the cylinder on which all C atoms are sitting (see Figure [TC]): 
note that implies that the resultant of this force system is null. 

As to couples, we guess that angle changes are equal at points C and F and 
at points D and E: 

Aa(C) = Aa(F), A/3(C) = A/3(i^); 

Aa(D) = Aa(£:), A/3(L>) = A/3(i;). ^ ^ 
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Figure 16: A double ABU: top-down view. 



and that, moreover, angle changes are opposite at points C and D and at points 
E and F: 

Aa{C) ^ ~Aa{D), A/3(C) = -A/?(i?); 

Aa{E) = -Aa{F), Ap{E) = -A^(F). ^ ' 

Combining pop and (|3ip . we arrive at the following set of symmetry-based 
assumptions: 

Aa{C) = Aa[F) - -Aa(i?) = -Aa(i?), 

A/3(C) = A/?(F) = -A/3(£;) = -A^p). ^ 

The equilibrium of our stick complex is guaranteed if the resultant moment with 
respect to any chosen pole of the forces and couples applied at points C, D, E, F 
is null and if, moreover, both the force and couple system acting on stick AC 
and the force and couple system acting on stick AD have null resultant moment 
with respect to point A. 

The moment balance for the double ABU yields the scalar condition: 

(ol cos \ (y, 
kg sin — — ; — — j A/3(C) — —{s sin — cosz9)/i . 

(33) 

The moment balances for the AC and AD sticks furnish a symmetry condition 
that we had not been able to guess, namely, 

Aa{A) = 0, (34) 
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and three scalar conditions: and 



ks £^ (A^(C) + A/3"(A)) = - (5 sin(7 -i)))h, 

( ct . \ /, . X A , • ctsin27^^,^, / . a . „\ , 

s cos — sm7 w — [ks sm7jAQ!(Cj + fc^ sm — ^ — — Ap(G) = — ssm — sin?y ft , 

V 2 / 2 smp \ 2 / 

fs cos ^ cos 7) w - (/c, cos7)Aa(C) +fc^ sin ^^^^^?^A;3(C) + fc^^^A/S^lA) 

V 2 / 2 sm p sm p 



sin ^ cos i9 ) ft. , 



2 



(35) 

where A/3"(^) denotes the change in angle between the AC and AB sticks. The 
developments necessary to arrive at these equations can be found in Section [8.31 
The solution of the system of ([33]) and (|35|) is: 



. a sm(t9 — 7) , , . , ^-1 OL Pa ■ „ 

V = — sm 4/1, with 4 :=cos — COS7 ^smw, 

2 sm7 2 s 

AQ;(C) = sm- cos-cosi?-— ^ ] ^hs ks, , , 

^ ^ 2^2 s sin7 y ^ ' " (36) 

A/3(C) = 0, 

A/3" (A) = sin/g '^"'^'^"^^ /is/fc,. 

sm7 



Note, in particular, ([SBjlg, another difficult-to-guess symmetry condition. Note 
also that it follows from the definition in (pG]) . with the use of ([2]), (|4]), and 
(|116p . that ^ = —2. Thus, the first two of ([36]) can be written in the following 
simpler form: 

. a sin(?9 - 7) , 
-y = 2 sm /i, 

''''^ (37) 
A ^ • sin(??-7) +sin7? , 

Aa(C) = 2smQ: ^ /is/fc^. 

sin 7 

It follows from the first relations in , (|115p , and , that the axial force 
in the sticks AC and AD is: 



/c • c =: iV(AC) = - AT (AD) = - f cos - cos(z? - 7) - 2 sin^ - ^^^^^ 

\ 2 2 sin 7 

thus, their axial deformation is: 

mAC) = -Ab(AD) = -(cos- cos('d - 7) - 2 sin^ - h/ka. 

\ 2 2 sin 7 / 

(38) 

Moreover, the axial force in stick AB turns out to be null, so therefore 

Aa = 0. 
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4.4 Radial loading of an ACNT 



Figure[T7]features a top view of a double ABU subject to radial forces of common 
intensity p, a load condition that we later assimilate to the radial pressure 
problem for a cylindrical shell. The front view is the same as in Figure [151 We 




1? \ 



C=D 




Figure 17: A double ABU subject to uniform radial forces: top-down view. 

presume the following symmetries, for the applied forces consequential to our 
Euler cut: 

/c = /d , fE = fp, 

fc ■ Ci= - fE ■ Ci , 
fc ■ C2= fE ■ C2, 
fcC3=fE- C3 = 0, 

and for the angle changes: 



(39) 



Aa{A) = Aa{C) =: Aa, A/3(A) = Af3{C) =: A/3, 

where Aa, A/3 satisfy 

Given and given that 

fA+fB = -2sin(ppc2, 

the resultant of the force system is null if 

-psmLpc2 +fc+fE = 0, 

i.e., if, in addition to the second of (HH]), 

sin if 



(40) 



fc ■ C2 



■P- 



(41) 



Only one more balance equation is needed, and we choose it to express the 
rotation equilibrium of stick AC with respect to point A: 



PC X /( 



C ^ 'a 



0, 



(42) 
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where, with the use of (|i(7|) . 



Tq = Tq c X d, Ta = kgAa |c x <i| ^ , 
Tfs = -T/s b X c, Tp = ksAl3 \c x 
fa = fa{-c) X {-d), fa ^ ksAa\c X d\~'^, 
= -T/3 e X (-c), fp = ksAl3\c X e\^^. 



(43) 



The algebraic manipulations in Section (|8.4I) allow us to conclude that 

sin a sin 



Aa = -- 



2 cos § sin 7 + (5 (a, /3, 7) sin a sin sin 27^ 
where, we recall, 



ps/fcs, (44) 



■^A Sill ^ 

S {a,l3,-f) = 7-T— 5 cos 7; 



2sin/3 



as to the axial deformation, we obtain: 



and 



Aa^ (S (a, 13, 7) sin (p)p/ka, 



Ab = 2 ( '''^^ 2" (^i'^T + T 7) )p/ka ■ 



(45) 
(46) 



4.5 Radial loading of a ZCNT 



Consider a double ZBU, subject to applied radial forces of common intensity p 
(Figure fT8|). and imagine it isolated from the CNT it belongs to by means of an 
Euler cut. 



P, 





B=F 



Figure 18: A double ZBU subject to uniform radial forces: top-down and front 
views. 
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We see that: 

fc ^fA^P VA, with rA ■■= - [sin cos - C2 j (47) 

(with reference to FigurelTSl-left. note that 4> — (tt — 7)/2). The force symmetries 
we presume are: 

fo = fs, fs = fp, fo ■ rA = fs ■ rA, fD-rA= -/b • rA, 

where := x C3; in addition, we presume that 

/d • C3 = /s • C3 = . (48) 

Under these presumptions, only two balance equations are needed for equilib- 
rium: the resultant of all forces applied at points A, B, C, D, E, and D and the 
resultant moment with respect to A of all forces and couples acting on stick AB 
must both be null (see Subsection 14.51 of the Appendix) . 

In view of (|47p . (|4.5p . and p8)) . the force balance equation reads: 



fB+fD+prA = 0, 
and is equivalent to its only scalar consequence: 

ifB+fD)-rA+P = 0, (49) 

or rather, 

sin -fs ■ ci + cos -fs ■ C2 = - p. (50) 
The moment balance for stick AB reads: 

PB ^ fB+ Ta +Tp + fo,+Ti3 = 0, (51) 

where, 

Ta^Tab X C3, Ta = ksAa\b X 03]"^ 

Tfi = -Tfj d X b, Tfj = ksAl3\d X b\-\ 

, (52) 

Ta^TabxC^, = fcs Aa | ft X C3 | , 

T/3 = —ffj m X {—b), fp — fcsA/3|d x b\^^ . 
The developments leading to establishing that 

^ \ P£ 

4 sin i (^2 - ^-4^^ (sin a cos 7 + sin a)^ ^^^^ 

AP ^f{a,l3,-f)Aa, /(a, /3, 7) = --^ (1 + C0S7) , 

2 smfj 

Aa = 0, 



and 



sma n 



are found in Section |83 
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4.6 Energetics 

The elastic energy stored in a deformation can be evaluated in a number of 
alternative ways. We find it convenient to do it with reference to bond units. 

4.6.1 The energy stored in A- and Z- CNTs under axial traction 

With reference to Figure [TUl the total energy stored in an ABU turns out to be: 



= fc„(A&)2 + fc,((Aa)2+2(A/3)2), 

(55) 

where Ab is given by ^3^^ ^-^d Aa, AjS by for a ZBU (Figure[TT]), we find: 
C^Icf" = \ka {{Aaf + 2{Abf) + \ks (2{Aaf + {Afif + 4 x ]^{Aaf + 2 x \[Apf 
= ha [{Aa)' + 2(A6)2) + ks (2{Aa)' + (A/?)^) , 

(56) 

with Aa, A6 given by d^T]) and Aa, A/3 by There are 4n ABUs and 2n 

ZBUs in a module, so that the energy stored per module is: 

^if=4nC/if^ [/If =2n[/|f^. (57) 

4.6.2 The energy stored in a twisted armchair CNT 

In this instance, the energy stored in an ABU is: 



1 



Aa)' + 2(A&(C))' + 2{Ab{D)^ + -k^ \2{Aaf + 2(A/3"(A))' 



TjABU _}_,, 

+ 2{AP\A)Y + {Aa{C)Y + {Aa{D)Y + (A/?(C))' + {Ap{D)) 

= 2ka{Ab{AC)f + k, [2 (A/3" (A))' + (Aa(C)) 

(58) 

where A;3"(^) is given by Aa{C) by and Ab{AC) by 

4.6.3 The energy stored in A- and Z- CNTs subject to radial loading 

The energy stored in an ABU is: 

= \ka{{Aaf + 2(A6)') + fc.(l + 2(/(a, /3,7))') (Aa)', (59) 

where Aa, A6, and Aa, are given, respectively, by (|45l). p6)) . and p4l) . For a 
ZBU, the stored energy is: 

C/|f^ = ka{Abf + (2+ (?''(a,/3,7))')fc.(Aa)', (60) 
where Ab and Aa are given by, respectively, (|54|) and (|53p . 
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4.7 Slenderness and thinness 

We here discuss two important properties of a SWCNT, its slenderness and its 
thinness, within the context of DSM. 

For the number of C atoms at our disposal (say, the largest number of 
atoms we can treat in a numerical simulation), those atoms can be arranged in 
a nanotube consisting of 'mt{N/NM) modules, where int(-) denotes the integer- 
part mapping. For Im the length of a module, we have the following formula 
for the length L of the nanotube in question: 



Consequently, that nanotube's aspect ratio - that is, its diameter /length ratio 
- is: 



where use has been made of ([2]) and pO)) . Therefore, a nanotube is slender if, 
roughly speaking, n? ^ N . Consistently, when matching theoretical predictions 
with computational results, convenient simulations should concern arrays of 
some N — 10^ v? atoms or more, as a glance to (|5Tt indicates. 

The chirality index n has to do not only with the slenderness of A- and Z- 
CNTs, but also with their thinness. A purely geometrical notion of thinness is 
easy to state within the framework of CSM: given a cylindrical shell of diameter 
2(Pq + e) and length 21, we term it thin if e/ <^ 1, that is to say, if its 
thickncss-to-diameter ratio is smallH Thickness, however, is an ill-defined notion 
for CNTs, especially when they are single-wall. This is why, in our opinion, a 
purely geometrical discrete notion of thinness makes no sense for CNTs. Be it as 
it may, on recalling relation ()22[) between radius and chirality, we can estimate 
as follows the 'geometrical thinness' of a (n,TO)— SWNT of arbitrary chirality 
and radius Pq: 



Thus, for an ACNT to be thin, the chirality index n must be larger than, say, 
10 or more, a rather rare occurence. We shall see in Section 16.31 that a more 
appropriate, not purely geometrical, notion of thinness does yield a decreasing, 
but less pronounced, dependence on n: shells mimicking the mechanical behavior 
of SWCNTs are generally not thin! 

*See I22| for a discussion of the non- geometric ingredients that a notion of thinness may 
contain. 





(61) 




(62) 
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5 Nanoscopic CSM. A SWCNT-Oriented Shell 
Theory 



In the majority of papers where CNTs are regarded macroscopically as shells, 
textbook theories induced from classic three-dimensional isotropic elasticity are 
used; consequently, the shell response is characterized in terms of two elastic 
modulilfl On the basis of the DSM modeling at the nanoscale of armchair and 
zigzag CNTs we here use, we think it better to view all CNTs, be they single- or 
multi-wall and whatever their chirality, as orthotropic cylindrical shells whose 
midsurface has a tangent plane coinciding with the orthotropy plane, so that the 
shell geometry agrees point-wise with the geometry intrinsic to the chosen type 
of material response (see Fig. [TOl where the three little cylinders suggest what 



Figure 19: The geometrical elements relevant to describe the material response 
of CNTs modeled as shells. 

probes one should cut out of the shell body in order to determine its material 
moduli). Such a shell theory has been fully developed in [S]; we here employ 
the simplest version fitting armchair and zigzag CNTs. 

5.1 Displacement and strains 

With reference to Fig. [20l given that (xi,!?, C) & (-^+0 x (0, 27r) x {-e,+e), 
the typical shell we consider as length 21, wall thickness 2e, and aspect ratio 
Po/l', we call it slender if p^/l <C 1. In the theory we import from |0, displace- 
ment fields inducing thickness changes and thickness shearing are excluded a 
priori. In addition, for our present purposes, consideration of axisymmetric so- 
lutions to axisymmetric boundary value problems suffices. These solutions are 

^There are, however, noticeable exceptions to this oversimplifying practice, such as |23| 
and [6]. 
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Figure 20: A portion of the model surface 5 of a right cyhndrical sheh. 



all independent of the circumferential coordinate and have the following form 
in physical components u<i>: 

u<i> — a<i> — (w' , u<2> = I 1 + — ) a<2>, u<3> = w, (63) 

V PoJ 

where a prime denotes differentiation with respect to xi, the only space variable 
from which all of the parameter fields a<i>,a<2>, and w, may depend. With 
the use of the linear strain measure 

E{u) = ^{\/u + \/u^), (64) 

it is not difficult to see that 

E{u)n = 0, 

i.e., that a shell of the type under study cannot change its thickness and that 
its fibers orthogonal to the referential midsurface S must remain orthogonal to 
it after any admissible deformation. The not identically null strain components 
are: 

E<\\> = a<i>' - Qw" , 

= K' ^ (65) 



E<22> = ( ( 1 + — ) ) W 



(needless to say, E<i2> — E<2i>). On defining the cross-section strain measure: 



E{n)--=:i r + E(n.(), (68) 



we have from (1651) that 



£^<ii> = a<i>', E<i2> = E<2i> = ^a<2>' , E<22> = p^^w . (67) 
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5.2 Constitutive parameters 

The mechanical response of the material comprising the cylindrical shells can 
be expressed in terms of five elastic moduli: Ei, E2,vi2,i^2i,G; the first two 
are Young-like moduli, the second two being Poisson-like, and the fifth being a 
shear modulus associated with the (ei, 62) pair of directions; the first four are 
not independent, because it must be that 

5.3 Axisymmetric equilibria 

Here below we record the equilibrium solutions of the axial traction, torsion, 
and pressure, problems; once again, the reader is referred to [S] for details. 

5.3.1 Axial Traction 

Suppose the only loading is a distribution of end tractions equivalent to two 
mutually balancing axial forces of magnitude 

F = (27rp„)p, with p = 0(e). (69) 

One finds that a<2> = and that, for a slender shell. 



w{xi) = -p„ 1/12 ^(e/Po) 
a<i>{xi) = (1 - 1^121^21 (1 - S{e/p^)))-^xi, 



(70) 



where 



H^/Po) — — 1+^ ; (71) 

^ log - :/l2i/21 



note for later use that 

1 1 + - 

lim — log f- =1, whence lim S{e/pg) — l. 

£/p„->-o 2— 1 ^ ^ e/Po-*-0 

5.3.2 Torsion 

When the only applied load is a distribution of end tractions statically equivalent 
to two mutually balancing torques of magnitude 

T = {2tt pl)t, with t = 0(e), (72) 

the axial and radial displacements a<i> and w vanish identically, whereas 

1 £~^t 

a<2>{xi) ^ -—^—xi. (73) 

p? 

As shown in [5], these results hold for whatever the value of the aspect ratio. 
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5.3.3 Radial pressure 

When the slender shell under study is subject to a uniform radial pressure 
w = 0{e), all the other applied loads being null, one again finds that a<2> = 0. 
Moreover, for a slender shell 

w[xi) = PoO[e/pJ—— and a<i> = -V2iO{e/pJ —— — xi. (74) 



6 The Effective Parameters of CNT-like Shells 

No matter whether we model a SWCNT as a shell or as a stick-and-spring 
complex, we may regard it as a cylindrical probe. Two parameters characterize 
the mechanical response of a probe in an axial traction experiment, and one in 
torsion; their verbal definitions are, respectively, 

axial load radial deformation ^^^^ 

axial deformation' axial deformation ' 

and 

axial torque , , 

STo = . . , ; (76) 

axial twist 

sax and Sto are stiffness measures, vax is a Poisson-like modulus of transverse 
contraction. Two additional parameters are needed to characterize the response 
of a SWCNT to radial loading, namely, the specific stiffness and Poisson-like 
moduli 

radial load axial deformation 



radial deformation radial deformation 

In this section we derive mathematical expressions for these five parameters 
within the frameworks of both continuous and discrete structure mechanics. 
The CSM expressions involve Ei,E2, 1^12, 1^21, and G, the five elastic moduli 
of the shell theory we adopt, as well as two geometric parameters, e and p^, 
the shell's thickness and model-surface radius. The DSM expressions depend 
on the nanoscopic spring-stiffness moduli ka and fcg, as well as on s, the C-C 
bond length, and Pq, the radius of the cylinder on which the C atoms arrange 
themselves (both for A- and Z- CNTs, Pq oc ns, recall ([2]) and ([TO))). We aim 
to evaluate the constitutive and geometric parameters of our nanoscopically 
informed theory of CNT-like shells solely in terms of fc^, and n. Our simple 
method consists in equating as many as needed of the continuous and discrete 
expressions we derive in the next two subsections for the various stiffnesses and 
contraction moduli listed above. 



^See for a full treatment of the traction and pressure problems for shells of arbitrary 
aspect ratio. 



27 



6.1 'Continuous' parameters 

When we regard a slender cylindrical shell as a probe subject to an axial load 
F, we set: 

41:=^ and - (77) 

E<ii> E<ii> 

for the effective axial stiffness and axial contraction modulus. With the use of 
(l67l) and dTOl), we then find: 



Ac) \ p A,. (c) _ 6{£/po) ,^ 

(78) 

where A(e) := ^^Trp^e is the area of the shell's cross-section. 
For a twisted shell of arbitrary slenderness, we set: 

e:=p-ia<2>', (79) 

where Q denotes the change in cross-section rotation angle per unit length 
caused by the application of the axial torque T. Accordingly, with the use 
of (|73p . we find for the effective torsional stiffness: 



GJ{s), Jie):=4TTple(l + ^) \ (80) 



where J{e) is the polar inertia moment of the cross section^ 

In case of a uniform radial pressure ro, we define for our slender shell-like 
probe an effective radial strain p^^w, in terms of which we set 



and .i^) 



-1 '^^^"^ -1 
Po W Po W 



for, respectively, the effective radial stiffness and the effective radial contraction 
modulus. With these definitions, relations ([74|) yield: 

ii1 = -2i. (81) 



(c) _ '2E2 ej Po (c) 



Remark. The dimensions of the stiffness parameters are: 

dim(s^],) = force, dim(s^]) — force x (length)^, dim(s^^) — force x (length)^^. 



^Recall that I I73II . and hence I I80II , holds whatever the slenderness of the shell under con- 
sideration. 
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6.2 'Discrete' parameters 

On resuming the spring-and-stick models of an A- and a Z- CNT assembled in 
Section 3, we set: 

(rf) _ F (rf) _ AF^yP^^ 

for the discrete axial stiffness and axial contraction modulus, where H'^^,P'^' 
and AH^ , AP^ are the undeformed height and perimeter of a module and their 
changes under the applied loads. While the first of these definitions parallels the 
first of (1771) in a self-explanatory manner, the second one requires a motivation: 
ours is that, for n sufficiently large and axisymmetric loadings, 

For an explicit evaluation of and v^fx-i '^^ proceed case by case. We find 
(Appendix, Subsection 18. 6^ : 

• (armchair case) 

Ad.A) 
■^Ax 



An 

kas (83) 



and 



1 1 /l I 3 A ka£. 

{d,A) _ 2tan^/jy fc. 

Ax ~ , \ -1 o 

o I 1 ( 1 I 3 \ k„s^ 



4 \ ' 2 tan2 (i 



(84) 



(zigzag case) 



and 



3 I 1 fca „2 



Ks (85) 



V 



1 , x/3(l+cos/3) fcgs'' 
(d,Z) _ 4 ' 8(5+cos^) ks 



(86) 



An alternative way to evaluate discrete strains and stiffnesses is to make 
use of Lame's theorem of work and energy ( 9 , Section 28). As is well known, 
this theorems guarantees that, in a linearly elastic equilibrium problem, the load 
work is twice the stored energy. Accordingly, for axial traction, 

FAH^^ = 2 U^i, with U^i = energy stored in a module. 

Hence, on applying definition (15^ ^ , wc find: 
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in particular, when is evaluated with the use of ([55]) and ([57|. follows. 

To estimate torsion stiffness, we invoke Lame's theorem and write, for a 
module, 

where T is the applied torque and Q is the twist angle per unit of length; on 
inserting this relation in (|80)) . we find: 

(d) _ J B 
^To — 9 TjM ■ 

For an ACNT, T = 2n hpQ, where h has been defined in Section 1431 (see Fig. [15] 
and equations On using ([58]), ^ and dST]), we find: 

^id,A) ^ 9(V3/2^^)n3 

cos(7 - ^) + 22iH!^ + 2 (2 -"-^^^^-'^ + 3 (i cos ^ - ^ t?^) ^ 

(87) 

In the case of radial applied forces, the Lame theorem yields: 

N''pwm = 2U^'^, (88) 

where WRd = radial displacement. It is the matter of a simple computation 
(Appendix, Subsection l6.2p to arrive at the following expressions for the discrete 
radial stiffnesses of a A- and Z- CNT: 

I J t\ 62 . 'i ' n 

4f = 7 V ^-^^ ' 

V3^(5 + 4cos27+^^^^j 

lMl7-8cos2/3+12cos7 + 3cos27)i^ , _i 

likewise, it is not difficult to show that 

^id,A) ^ ^ 2^-ll + 2cos2/3(4-^)-3cos27 

sin2/3 4(3^+4) + 2cos7(64 + 3cos7ii|¥p' ' 

and 

^k'/^=-^ + ^ ^ TTTY. (92) 



311 ■ ^ 



(6 + 9^) 
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6.3 Elastic moduli, thickness, and radius, of CNT-like 
shells 

i. Contraction moduli . We choose the cylinder axis of a CNT-like shell parallel 
to the roll-up axis of an ACNT, and require that 



Rd 



Rd ' 



(93) 



whence, by ((81| g and (|9T|1 . we obtain that 



1^21 



3 2^ - ll + 2cos2/3(4- ^) - 3cos27 



sin /3 4 



(3^ 



4) + 2cos7(64 + 3cos7ii|^ 



(94) 



Likewise, on choosing the cylinder axis of a CNT-like shell made of the same 
orthotropic material parallel to the roll-up axis of a ZCNT, we get: 



(d,z) 
^12 = '^m 



3 1 



(95) 



The dependence of these two Poisson-like moduli on the size parameter n is 
shown in Figure [HJ note that 



0.35 

!^12 ••• 0.3 
^21 OOO 

0.25 

0.228 
0.2 
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Figure 21: The contraction moduli 1^21 (o) and 1^12 (•) of ACNT-like shells. 



lim V21 = lim U12 =: = 0.228. (96) 

ii. Radius . With a view toward determining the effective radii and of A- 
and Z- CNT-like shells, we observe that, by the use of ^S^li, (|HIl)i and (ffSjl o. 
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and for 



to hold, we find: 



El 



2^ Ax 



El 



'^12 E2 



(c) 
V2I 



whence, thanks to (|8iP o. we arrive at 



(c) (c) 
^Rd ^Rd 



^ Po = 



\ 



^ '^Ax ^Rd 

27r 

^Rd ^Ax 



We then set: 



\ 



1 s 



{d,A) 
'^Rd 



27r J'i.'4) , ' 



'fid 



1 S 



{d,Z) {d,Z) 



Ax 



Rd 



Xd.z) (d.z) 

'Rd '^Ax 



(97) 



Botlr effective radii depend on the size parameter n in a complex manner. But, as 
Fig. l22l shows. neither of them differs much, whatever n, from the corresponding 



Ya oc 



1.1- 



1.05- 



1 - 



0.95 



0.9 
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Figure 22: Ratio of effective and geometrically necessary radii of A-(o) and Z-(») 
CNT-like shells. 

value of the geometrically necessary radius; moreover, as n grows bigger they 
both tend to it: 

lim Po/Po = lim Po/Po = 1- 
iii. Thickness . We begin by laying down the following equality: 

H^^/Po) ., _,.{d,A) 



1- 1^121^21 il-Sie^/p^)) 



■ ^12 = ''Ax 



(98) 
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where the values for 1^21, J^i2 and are chosen as agreed in i. and iL, and give 
it the role of an implicit equation that we solve for the unknown 

As to , we proceed as follows. Firstly, we equate the continuous and 
discrete radial stiffnesses of both A- and Z- CNTs: 



2^26 -Sj^^ p^d(£ Ip^), 
Zhie - s^^ p^d{e j p^) 

(cf. (EUi, ([Ml), and (inni)), whence, with the use of ([55]), we arrive at 



(99) 



,, M^^') si( ^Z I ^Z\ 
£ _ ^21 Sflrf Po °\'^ IPo) 



-A 



yi2 4^/) P^Ke^lP^)' 



Finally, given that (|98| implies that 



' Ax 



1 - 1^121^21 



'5(^"/^o) = ^- (100) 

and that a completely similar relation holds for (5(e'^/pf ), we obtain: 

(yd.Z) {d,Z) , z 

, ^A ^Rd ^Ax 1 ~ ^^21^^^^ Po , . 

(d,A) (d,z) . (d,z) pA ■ y^^^> 

^Rd ^Ax 1 - ^^■^^Ax ^° 

Perusal of ([98]) and (jlOip yields for the effective thicknesses and the 
dependences on the chirality index n shown in fig. [23l 

iv. Extension moduli . Once effective radii and thicknesses have been computed, 
one way (there are others) to find the constitutive moduli E\ , E2 is to combine 

(c,A) _ U.A) , ^(c,A) _ ^(d,A) 



78]) and ()100p . under the assumption that sj^'^ = s)^^ and Vj^^ ~ Ax 



so as to obtain: 



a germane relation holds for E^, namely, 

E.i.2e^)^ ''fl (103) 

The dependence of both i?i(2£"*) and E^i^e^) on n is shown in Figure 
Figure YIM shows how the values of both E\ and E2 decay with n, their difference 
becoming smaller and smaller. 



*With the use of I I71II . it is the matter of a short computation to give I I98II the following 
form: 

1 1 + a: _ 1^21 
for X the unique positive solution of this equation, wc have that e-^ = x. 
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Figure 23: The effective thickness of A-(o) and Z-(») CNT-like shells. 



V. Shear modulus . By equating the continuum and discrete evaluations (j80p 
and (|57| of the torsional stiffness, we find: 

(d) (d) 

(see Figure !^ . 

vi. Consistency checks. The correctness of our computations can be checked in 
a number of ways. For one, (|68p. a key feature of our description of the consti- 
tutive response of orthotropic shell, is satisfied exactly, to all practical puposes, 
by the quadruplet (i?i, i?2, J^i2, ^^21) computed in the preceding subsections by 
equating continuous and discrete stiffnesses and contraction moduli. In fact, by 
an appropriate re- use of the computations made to construct Figures [2l1 and [24l 
one can show that, whatever the index n, E1/E2 — 1^12/1^21 ~ 5 x 10~^^. 

For two, the thickness ratio je^ can be evaluated in a way alternative to 
the one implicit in (|23p. namely. 



^~^A„Z ^Rd ^Ax i-~V2lVAx m (^nr,\ 



'Aa; 



^^21 1 - V12V\ 



with numerically identical results. Interestingly, no matter whether (j23p or (llOSp 
is used, one has that 

lim -4=1, 

that is to say, the thickness of A- and Z- CNTs of larger and larger radius tends 
to be the same. 



{ToUl is arrived at without presuming that holds. 
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Figure 24: The extension moduli of CNT-like shells. 



7 Conclusions 

Wc have presented a proeedure to determine the constitutive moduli and the 
geometric parameters to be associated with a single-wall carbon nanotube of 
given chirality in order to model it as a linearly elastic orthotropic cylindrical 
shell; our findings are summarized in the following table, for various values of 
the chirality index n, where we have set h-^ := 2 e"^, := 2s^ . 



n 


6 


10 


20 


El (TPa) 


0.832 


0.736 


0.709 


E2 (TPa) 


0.784 


0.729 


0.706 


1^12 


0.260 


0.241 


0.232 


1^21 


0.242 


0.233 


0.230 


G (TPa) 


0.424 


0.359 


0.307 


(nm) 


0.390 


0.446 


0.466 


Po (nm) 


0.402 


0.675 


1.354 



n 


6 


10 


20 


El (TPa) 


0.784 


0.729 


0.706 


E2 (TPa) 


0.832 


0.736 


0.709 


1^12 


0.242 


0.233 


0.230 


V21 


0.260 


0.241 


0.232 


G (TPa) 


0.424 


0.359 


0.307 


(nm) 


0.388 


0.446 


0.466 


Po (nm) 


0.245 


0.399 


0.787 



Table 1. Elastic moduli, thickness and radius of CNT-like shells (A-, left; Z-, right). 

As anticipated in our introductory section, the list of constitutive and geo- 
metric parameters appropriate for a SWCNT of arbitrary chirality pair (n, m) 
can be deduced for either one of the two lists in Table 1. A feature of our pro- 
cedure is that each of these parameters, in addition to chirality, depends only 
on the two material constants modulating the harmonic approximations of the 



35 



TPa 



G»»» 



1.5 



0.5 - 



10 



Figure 25: The shear modulus of ACNT-hke shells. 



nanoscopic bonding energies and on the length of a C-C bonds 

Here is an interesting conclusion that can be drawn from our results in the 
so-called 'graphene limit', that is, for index n — >■ oo. We have that 



E1/E2 



4'/^ e^ l-i^i2i^2i{l-6{e^/p^)) 



1, ^ 1, h^/h' 



moreover, 



lim Gh"^ =: G°°h°° = 0.143 TPa x nm, 

n— J-OO 

lim Eih"^ =: E^h°° = 0.332 TPa x nm. 



= 0.228 



(recall ([M)) ). Now, it is to check that 



2(1 



= 0.135 TPa X nm^ G°°h° 



whereas equality should hold when graphene, modeled within linear plane elas- 
ticity, is regarded as isotropic. 

Another feature of our study is that it allows for determining in a mutually 
consistent manner the values of both wall thickness and Young modulus of a 
SWCNT. What value to take for thickness has been an issue in the nanotube 
community for long. Needless to say, an evaluation of effective thickness is 
necessary to apply any shell theory, ours being no exception. In the following 



^''These three nanoscopic bits of information can be put together under form of a dimen- 



sionless characteristic number, namely, — ^ 
stretching and angle-changing energies. 



13.101, measuring the relative importance of 
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table, a subset of Table 1 of [TT], values for wall thickness and Young modulus 
of a SWCNT are listed, as determined by a number of different procedures; we 
have added a last column, for the values of the tensile rigidity [TT], that is, the 
product Young's modulus x wall thickness. 









Wall 


Young's 


Tensile 


Authors 




Method 


thickness h 


modulus E 


rigidity Eh 








(nm) 


(TPa) 


(TPaxnm) 


Lu 1171 




Molecular dynamics 


0.34 


0.974 


0.331 


Hernandez et al. |10| 


Tight-binding molecular dynamics 


0.34 


1.24 


0.422 


Li and Chou |15l 


Structural mechanics: stiffness 


0.34 


1.01 


0.343 






matrix method 








Jin and Yuan |13| 


Molecular dynamics 


0.34 


1.238 


0.421 


Yakobson |30l et al. 


Molecular dynamics 


0.066 


5.5 


0.363 


Zhou et al. [32 




Tight-binding model 


0.074 


5.1 


0.377 


Kudin et al. 




Ab initio computations 


0.089 


3.859 


0.343 


Tu and Ou-Yang 


26] 


Local density 


0.075 


4.7 


0.352 






approximation model 








Vodcnitcharova and 


Ring theory continuum mechanics 


0.0617 


4.88 


0.301 


Zhang [271 












Pantano et al. |21| 


Continuum shell modeling 


0.075 


4.84 


0.363 


Wang ct al. [29] 


Ab initio calculation 


0.0665 


5.07 


0.337 



Interestingly, while the reported values for wall thickness and Young modulus 
display a large scattering - the so-called Yakobson's paradox [3S] - the tensile 
rigidity has a very small variance, irrespectively of the method used to evaluate 
the factors. We regard this evidence as an implicit resolution of the paradox, 
in that it makes evident what one can sensibly try and evaluate: the stiffness 
parameters defined in Section |6l or others of the same sort, in two ways to be 
juxtaposed, discrete and continuous. 

Note that the values furnished by our theory for the tensile rigidity comply 
with those in the table whatever the chirality index; Figure [26] summarizes our 
findings. 

Finally, it is worth noticing that, in the light of our present theory, two 
notions of thinness can be given for a SWCNT. The first, geometric thinness, 
is expressed by formula ([5^ : it depends exclusively on the chirality index n, by 
way of the ratio of the length s of the C-C bond and the geometrically necessary 
radius pg. The second, effective thinness, is e/ p^; here both e and depend in a 
complicated way, in addition to s and n, on the constitutive nanoscopic stiffness 
parameters ka and kg. Figure [27] shows that the effective thinness, while being 
greater than the geometric one, approaches zero as n grows in the same n~^ — 
way as the latter; it also shows that for small radii CNT-like shells are far from 
being thin. 

8 Appendix. Ancillary Computations 
8.1 Computations relative to Section 14.11 

Note that 

|f) X c| = sin ^( = |c X e| , \cxd\=sina, 
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Figure 26: The tensile rigidity of A-(o) and Z-(») CNT-like shells. 



two relations that will be useful to derive l |106| l here below. Combination with {TtJ-j^, | |18^ -I |21| |. and 
^221 , allows us to give l\20l the following form: 



(s/c+2r„d + r^(e-6)) = 0: 



(106) 



thus, this equation has only two scalar counterparts. Now, with the use of ^5^^^ J16| l. l |17t . and 
((22j, one obtains: 



moreover 



< fc ^ ifc ■ C3) COS — sin 7C1 + y{fc ■ C3) cos — cos 7 - - (/b ■ ci) sin — ) C2 + 

1 / ^ a 

+ -(/b ■ ci) cos — sm7C3; 

it follows from l[5j and {tJ that 



(107) 



c X d — — sin a. (sin 7 ci + cos 7 C2 ) , 

a a 
c X e — — sin — (sin 27 ci + cos 27 C2 ) — cos — sui 7 C3 , 

a a 
c X o — sin — C2 — cos — sin 7 C3. 



(108) 



Consequently, in view also of jlSl l. | |19| > and | |21| >. the vector equation (IIO6I 1 yields the following 
scalar system: 



_ , sm 

scos — (/c ■ C3) - 2A;s Aq - ^ — ^ cos 7 ^ 0. 



sin j3 



(109) 



Thus, as to applied forces, wc have: 



fs = 0, fc = -fo = fc3, 



(110) 



These developments lead to ]23t : ( 124l > is arrived at by combining the second of JlOOt with and 
ig and fTTot . 
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Figure 27: Geometric (•) and effective (o) thinness. 



8.2 Computations relative to Section 14.21 

On taking il2l into account, one finds that 



'^Pb X /b = -(/b • C2) cos QCi + {^fs ■ ci)/ cos Q + - (/c ' C3) sin a ) C2 
+ (/b • C2) sin QC3; 



moreover, 



b X d — — sin a cos a sin 7C1 — sin a cos a(l -f cos 7)c2 + sin ck sin 7C3 , 
6 X m — — sin a cos ct sin 7C1 — (1 + cos 7) sin a sin 7C2 + sin^ a sin 7C3 , 
6 X C3 — — sinaC2. 

With the use of these relations, equation ]26t yields the following scalar system: 

/b • C2 = 0, 

(/b ■ ci)s COS a H (fc '^3)5 sin a — 2fcs (1 + cos 7) sin a cos a — 2A:^, Aa — 0. 

2 sin /3 



(111) 



Now, the first of fTTll and the last of ( 1251 1 imply that fjj - a X C3 — 0; with this and the penultimate 
of l |25l l, the force balance | |16l l imply that fs-ci— 0. Summing up, the force system is: 

fs^fo^ -fc ^ -/C3, 
whence the result in l]27|t . Moreover, ijllU lo reduces to 

1 fs AP ^ 

sm Q — 2 ( 1 + cos •y) sm a cos a — 2Aa — 0, 

2 A;^ sin^^ ^ 

whence, with the use of wc obtain ( 1281 1 . 
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8.3 Computations relative to Section 14.31 

(i) Preparatory developments. With a view toward formulating the equihbrium conditions for the 
AC and AD sticks, wc supplement | |32t with the following further guess: 

A/3''CA) - -A^'*{A), (112) 

where A/3" {Aj3'^) denotes the change in angle between the AC (AD) and AB sticks. Moreover, 
this time in preparation for formulating the first moment balance, we pick the pole P shown in 
Figure fl5l and, with the help of Figure f 161 we specify as follows the relative position vectors of the 
points where forces are applied: 



with 



PC — h rc, Pd — —h -\- rc , 
Pe — h ve, Pf — —h + rE, 



h — s sin — c-i , 

a 

rc — —po (sin -dci ~\- cos ^9c2) — ta + s cos ~ 
rE — Pq (sin 'dci — cos t9c2) — r^+s cos — a , 



(113) 



(114) 



and 



-/3q(cos (pci + sin if 02); ru — Pq(cos i^ci — sin (pc2), 
- cos 7C1 + sin 7C2 , a — cos 7C1 + sin 7C2 , 



pQ — — n s. <p — arccos — , 
27T 3n 



After some manipulations, we find: 

27r Of 27r a 

sin 1? — cos if H cos — cos 7. cos 1? — sin ip cos — sin 7, 

3n 2 ' 3n 2 

whence 

tc — cos 'dci — sin ^3c2 , is — cos ^9ci + sin i9c2, (H^) 

with 

1? ^ i9(n) :^ arcsin f — (l + cos ) V (116) 

V 3n V 2n / / 

we note here for later use that sin('(9 — 7) > 0, for all n > 2. 

(ii) Moment balance for the double ABU. The resultant moment of the force and couple system 
with respect to pole P reads: 

mp + T^(C) + Tf>(C) + f „ (13) + [D) + T^(E) + Ti3{E) + T^{F) + T^(F) = . (117) 

Here, mp is the moment of the force system with respect to pole P: 

mp ^ PC X fc -\- Pd X fD ^ Pe x fs -\- Pf X /f , 

for which, with the use of ([Sj, ( |29t , and JllSt -l fllSt , we find: 

mp — 4(/i s sin — cos -\- v pQ sin 'd^ C2 ■ (118) 



Ml 



oreover, 



-r„(C) = 


r„(C) (-C) X 


m, f^(C) = ksAa{C)\cxm\ 




TtAC) = 


-f„{C)ex ( 


-c), T3(C) = A;^A/3(C)|e X c| 


-1 


r„(z3) = 


-T„(D)(-C) 


X m, Ta{D) = fc3Aa(_D)|c X 


m| , 


niD) = 


ff,{D){-d) > 


c e, T^(D) = fc,.A/3(£))|(i X e|" 


-1 


r^{E) = 


— T„ (iJ) m' X 


(-c'), ?„(£;)= fc,Aa(£;)|m' 


X c'r^ 


Tp{E) = 


f0{E){-c') ; 


< (e'), ffiiE) = k^Al3(E)\c' X 




T^(F) = 


f Q (F) m' X ( 


-c'), Tq(F) = fcsAQ(F)|m' X 






-T;3(F)(e') ; 


< (m'), ff)(F) = fcsA,3(F)|e' > 





(119) 
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where 

e — — cos 27 ci -f sin 27 C2 , e — cos 27 ci + sin 27 C2 ; 

, a / a / a / a , , (120) 

c — cos — a + sin — C3, rf — cos — a — sin — C3; m — —d, m — —d . 
2 2 2 2 

Making use of jllSl l , l |119l > and ( |120t . it is not difficult to see that, as a consequence of the symme- 
tries embodied in l[32j, jllTI ) is equivalent to the scalar equation {33j 

(iii) Moment balances for the AC and AD sticks. For A the common pole, we request that 

PC xfc+T^{C) + Ti3(C) + T^{A) + ri(A) = 0, 

_[ (121) 
PD X fD+T^{D) + Ti3{D) - T„(A) + T^(A) = 0, 

where the terms not listed in ( I119t arc defined as follows; 

T^{A) = T^(A) cx d, To(A) = ksAa{A)\c x dr\ 

T^(^) = -T^(A)h X c, r^(A) = fceA/3"(A)|6 X c|-\ (122) 
T^(A) = rjiA) dx b, rf^iA) = k,Af3''{A)\b x dl'V 

Then, equations J12H can be rewritten as 

cx (s/c + (r„(C)+r„(A))d-r3(C)e + T^(A)6) =0, 
dx i^sfo + {To^iD) - Ta.(A))c - ff,{D)e + rj{A)b) =0. 

Note that Aa(A) is the fifth and last unknown in our equilibrium problem, the other four being v, 
A.a{C), A/9(C), and A/3"; the determining equations are l |33t and the four scalar consequences of 
< |123l l : by subtraction of one of ( 1123t from the other, we find both ]34| l and ]35| l ; by summation, 
we find ( 1351 1 2 ..-^ ■ 

8.4 Computations relative to Section 14.41 

Combining l\42l and l |43| l. we have: 

c X (sfc + 2Tad + r^b — r^e) — 0, 

an equation that, on recalling the definitions of the vectors involved, can be given the following 
form: 

sfc + 2ra d Tf^b — Tf^e — Ac, with A — — 2rn , 

or rather, equivalently, 

sfc + 2r«(c + rf) + r^b- r^e - 0. (124) 
The scalar consequences are two: 

s{fc • Ci) — 4 cos — cos 7 Ta + + COS 27 — 0, 

a 

s{fc ■ C2) + 4 cos — sin 7 — sin 27 — 0. 
If, in the second, we make use of l |41| > and ([9j, we find ( |44| >; with this, the first gives: 

fQ - ci — — — cot 7 sin p . (125) 

As to axial forces and strains, in the a-typc stick AB they are, respectively, 

- 2fc ■ b ^: N(AB) ^ cot 7 sin ^p, (126) 
whence we obtain l |45| l: in all other b-type sticks they arc: 

f(j ■ c :— N{AC) — 2 ( ""^^ ~2 ^^^"'^ ^ "^^ ^' C-'-^'^) 

whence ( 146| l follows. 
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8.5 Computations relative to Section 14.51 

Combining | 151I > and ]52t , wc have: 

'> X (s/b + 2TqC3 + Tpd — ff)in) = 0, 

an equation equivalent to 

sfB + 2t„(c3 +2b) + Tuid - m - 2b) = 0, 
whose two scalar consequences arc: 



/ ^(/b ■ ci ) + 4:Ta sin Q — 2r^ (sin a cos 7 + sin a) — 0, 

\ /S • C2 = 0. 



(128) 



Equations l |50t and ( II28I I2 allow to conclude that 

1 

2 sm ^ 

With this, ifml -^. and {TUl, we obtain i53l . 

Finally, while it is obvious that Aa — 0, for A.b wc have that 



1 sin a p , , 

8.6 Computations relative to Section 16.21 

To arrive at relations l |83| l and I l84| >. one starts from 
F = 2nf 

H^"' =4bsin-. P^"' =2n(a + bcos- 
2 V 2 

^^AM ^ 4gjjj -A& + 2()cos - Aa, AP"^"' = 2n ( Aa + cos -Afc - ibsin -Aq 
2 2 V 2 2 2 

(here, a. b denote the common undeformcd length of sticks of type a and fc, respectively), whence 

(3 + - 1 + ^ ^ ' 



H^M 2^/3^ i\ 2tan2^y fc, J k^s 

(1- - 1-' ^ a \ J 



P^" 2V3V 4V 2tan2/3y fe^ ^fcas' 
Quite similarly, equations (I85t and i86t are arrived at on setting: 
F = n/ 

H —2{a — bcosaj. P — inosma, 

AH^" = 2(Aa - cosaAfc + fcsinaAa), AP^''' = 2n(sin aAfc + 6 cos aAa), 

whence 

AH^" _(Z 3(l + cos/3) fc„s^\ / 



if^J>-f 8(5 + cos/3) fcs / feaS ' 

AP^*^ _ / 1 73(1 + cos^i) /c„s^ \ / 



pz'^'-f \ 4 8(5 + cos /3) fc, / kaS 
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It follows from ]88t that 
where U^^^ is dcduciblc from |58j by using (I45t and ( I46t : 



A „ ABU 



5 + 4 cos 27 - 



"'^■^ 16/Ca sin-^ 7 ^ ' o I g /" cos 7 ^ ^ ks 



Thus, 



A P sin V ( , , . „ , 6 kas' 

"Rd = r^— 5 + 4 cos 27 H — - 

Sfc,. sm2 7 I 8 + 3(g^) 
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